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Abstract 

J  A  computer  program  was  written  in  the  MUMPS  language  to  design 
filters  for  use  in  cancer  radiotherapy.  The  filter  corrects  for 
patient  surface  irregularities  and  allows  homogeneous  dose 
distribution  with  depth  in  the  patient.  The  program  does  not  correct 
for  variations  in  the  density  of  the  patient. 

The  program  uses  data  available  from  the  software  in  Computerized 
Medical  Systems  Inc.'s  Radiation  Treatment  Planning  package.  External 
contours  of  General  Electric  CAT  scans  are  made  using  the  RTP 
software.  The  program  uses  the  data  from  these  external  contours  in 
designing  the  compensation  filters.  The  program  is  written  to  process 
from  3  to  31,  1cm  thick,  CAT  scan  slices. 

The  output  from  the  program  can  be  in  one  of  two  different  forms. 
The  first  option  will  drive  the  probe  of  a  CMS  Water  Phantom  in  three 
dimensions  as  if  it  were  the  bit  of  a  routing  machine.  Thus  a  routing 
machine  constructed  to  run  from  the  same  output  that  drives  the  Water 
Phantom  probe  would  produce  a  three  dimensional  filter  mold.  The 
second  option  is  a  listing  of  thicknesses  for  an  array  of  aluminum 
blocks  to  filter  the  radiation.  The  size  of  the  filter  array  is  10 

/  f  »  • 

-inches  by  10  inches.  The  Printronix  printer  provides  an  array  of 

•  r  '  * 

blocks  1/2  tnctrby  1/2  inch- with  the  thickness  in  millimeters  printed 
inside  each  block.  , 


COMPUTER  DESIGNED  COMPENSATION  FILTERS 


FOR  USE  IN  RADIATION  THERAPY 


I.  Introduction 


Surface  irregularities  of  cancer  patients  make  radiation  therapy 
more  difficult.  The  surface  irregularities  cause  the  treatment  beam 
to  pass  through  different  amounts  of  tissue  to  reach  the  depth  of  the 
tumor.  The  beam  is  attenuated  more  in  the  regions  where  it  has  passed 
through  more  tissue.  This  causes  a  skewing  of  the  isodose  curves 
which  cause  difficulty  in  treatment  planning  (Ref  1,2, 3, 4,6,8).  The 
patient  can  be  made  to  approach  the  ideal  of  a  flat  surface  orthogonal 
to  the  beam  by  the  introduction  of  a  filter  to  compensate  for  tissue 
deficits  at  the  patient's  surface. 

The  tissue  compensator  cannot  be  placed  directly  on  the  patient's 
surface  because  the  electrons  emitted  from  the  compensator  would  burn 
the  patient's  skin  during  treatment.  The  tissue  compensator  must  be 
at  least  15  cm  from  the  patient's  skin  to  prevent  burning  (Ref  3,4). 
The  displacement  of  the  tissue  compensator  above  the  patient  makes  it 
necessary  to  reduce  the  lateral  dimensions  of  the  compensation  filter 
to  correct  for  beam  divergence  (Ref  1,5,6).  See  Figure  1.  The 
vertical  dimension  of  the  filter  must  be  changed  by  the  ratio  of  the 
attenuations  of  the  compensator  material  and  the  tissue. 
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Figure  1.  This  is  an  example  of  the  concept  of  a  tissue 
compensator.  S  is  the  radiation  source  and  P  is  the  central  axis 
point.  C  is  the  filter  for  the  compensation  of  the  tissue  deficit 
D. 
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There  have  been  many  methods  of  designing  and  constructing 
compensation  filters  used  in  the  past.  Most  of  these  methods  were 
time  consuming  and  uncomfortable  for  the  patient.  One  such  method 
involved  the  use  of  the  "rod  box"  to  measure  the  topography  of  the 
patient.  This  involves  clamping  the  many  rods  coming  out  of  the  box 
into  position  so  that  the  rods  are  touching  the  patient's  surface. 

Then  the  comparative  lengths  of  the  rods  are  used  in  the  design  of  the 
filter.  Another  method  used  a  router  to  carve  out  a  styrofoam  filter 
mold  while  a  pointer  connected  to  the  router  was  run  over  the  surface 
of  the  patient  (Ref  1). 

The  idea  is  to  find  a  way  of  designing  and  constructing  a 
compensation  filter  which  is  less  time  consuming  and  is  more 
comfortable  for  the  patient  than  current  methods.  Cancer  patients 
routinely  undergo  computerized  axial  tomography  (CAT)  prior  to 
treatment.  These  CAT  scans  could  be  used  for  computer-designed 
compensation  filters  using  the  data  available  from  the  Computerized 
Medical  Services  Inc.'s  (CMS)  Radiation  Treatment  Planning  software. 

Problem 

The  problem  undertaken  in  this  study  was  to  use  the  CMS  computer 
facilities  available  at  the  Wright-Patterson  AFB  Medical  Center  to 
design  a  compensation  filter.  The  computer  would  use  the  external 
contours  of  the  CAT  scans  to  “*  -hip  Ihe  primary  aim  was  to  set  up  a 
MUMPS  language  program  to  drive  the  probe  of  the  water  phantom  as  if  it 


were  a  router  bit  being  used  to  carve  out  a  styrofoam  filter  mold  (see 
Appendix  A).  The  secondary  aim  was  to  have  a  hardcopy  printout  for  an 
aluminum  block  filter  to  be  used  until  the  routing  machine  could  be 
built  (see  Appendix  F).  The  program  would  also  need  the  versatility 
to  work  for  various  compensation  materials. 


Scope 

The  scope  involved  getting  the  program  to  run  for  a  twenty-one 
slice  sequence  of  CAT  scans.  The  CAT  scans  imaged  the  Rando  Phantom, 
a  human  skeleton  encased  in  tissue  equivalent  material  (see  Figure 
5).  Testing  of  program  output  was  limited  to  driving  the  water 
phantom  probe  in  what  appeared  to  be  the  proper  manner,  geometric 
evaluation  of  the  filter  and  film  dosimetry  using  the  compensation 
filter  with  treatment  of  the  Rando  Phantom  . 


Assumptions 

Two  major  assumptions  were  made  to  simplify  the  problem.  The 
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first  was  that  the  human  body  has  a  uniform  density  of  1  gm  /  cm 
and  the  second  was  that  effects  of  scattering  by  the  filter  and  the 
filter  holder  could  be  ignored. 

General  Approach 

The  general  approach  was  to  retrieve  the  position  coordinates  for 
the  external  contours  from  the  CAT  scan  slices  and  manipulate  this 
data  so  it  could  be  used  by  the  routing  machine  driver  routine.  Then 
the  data  was  further  manipulated  for  use  in  the  printout  of  the  two 
dimensional  array  that  would  be  used  to  design  the  aluminum  block 
filter. 
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Sequence  of  presentation 

A  more  detailed  analysis  of  the  problem  and  the  development  of  the 
program  is  discussed  in  the  next  section.  The  third  section  covers 
the  validation  of  the  project.  The  fourth  section  discusses 
conclusions  and  results.  Appendix  A  contains  information  concerning 
the  Water  Phantom.  Since  most  people  lack  familiarity  with  the  MUMPS 
language,  a  short  guide  is  in  Appendix  B.  Appendix  C  contains  the 
program  documentation.  Appendix  D  contains  the  program  itself. 
Appendix  E  contains  the  aluminum  block  filter  printout  . 


The  data  contained  in  the  stored  external  contours  of  the  CAT  scan 
slices  are  in  X  and  Y  coordinates  for  each  point  on  the  contour.  Each 
of  the  X  and  Y  values  is  an  average  over  the  one  centimeter  thickness 
of  each  CAT  scan  slice.  A  successive  series  of  these  CAT  scan  slices 
along  the  Z  axis  will  be  used  by  the  program  to  design  the  three 
dimensional  compensation  filter  (see  Figure  2). 

The  RTP  software  and  the  program  itself  must  be  copied  onto  the 
fixed  disc  of  the  disc  drive  unit.  The  Master  Patient  File  (MPF)  disc 
must  then  be  placed  on  the  disc  drive  unit.  The  General  Electric  CAT 
scan  slices  must  be  copied  from  the  magnetic  tape  onto  the  MPF  disc. 
The  CMS  software  allows  this  to  be  done  fairly  easily.  An  empty  MPF 
disc  has  room  for  only  twenty-four  CAT  scan  slices.  If  more  than 
twenty-three  CAT  scan  slices  are  to  be  used  in  the  filter  design  one 
must  copy  only  the  first  twenty-three  onto  the  MPF  disc,  then  external 
contours  are  made  of  these  slices  using  the  RTP  software.  These 
external  contours  must  be  initiated  with  the  light  pen  located  at  the 
bottom  of  the  CAT  scan  slice  (this  allows  a  later  simplification  of 
the  program).  After  these  contours  have  been  made  and  stored  with 
their  contour  descriptions,  the  CAT  scans  on  the  MPF  are  erased  so  the 
remaining  CAT  scans  can  now  be  copied  onto  the  disc  and  have  their 
external  contours  generated  and  stored  on  the  MPF. 


Figure  2.  The  figure  Is  an  example  of  the  treatment  area  of  the 
patient  as  a  composition  of  successive  CAT  scan  slices  and  shows 
the  coordinate  axes  at  the  patient. 
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The  X  and  Y  coordinates  of  the  central  axis  point  must  be 
determined.  The  central  axis  point  is  that  point  directly  below  the 
beam  source.  The  coordinates  can  be  found  with  the  existing  RTP 
software.  Using  the  External  Beam  Planning  option  a  treatment  beam  is 
positioned  so  that  It  passes  through  the  central  axis  point.  The 
hardcopy  printout  of  the  treatment  plan  gives  the  values  in 
millimeters  for  the  X  and  Y  coordinates  of  the  central  axis  point. 
These  values  are  used  later  in  the  program. 

The  program  must  be  divided  into  smaller  programs  because  of  the 
limited  size  of  the  buffer  in  the  CMS  computer  system.  Each  of  these 
smaller  programs  calls  up  the  next  program  in  the  sequence.  The  size 
of  the  buffer  also  limits  the  number  of  local  variables  which  can  be 
carried  along  from  one  linked  program  to  the  next.  This  requires  that 
most  variables  must  have  their  values  stored  in  global  variables  on 
the  disc  rather  than  in  the  buffer.  This  increases  the  amount  of  time 
necessary  to  run  the  program  . 

Program  Development 

The  first  problem  to  be  addressed  is  the  reading  of  the  values  for 
the  X  and  Y  coordinates,  for  each  point,  on  each  of  the  external 
contours  from  the  MPF  disc.  First  the  number  of  data  points  contained 
in  the  external  contour  must  be  read.  The  program  then  reads  the  X 
and  Y  coordinates  for  each  of  the  points.  The  coordinate  system 
origin  is  within  the  CAT  scan  slice  and  can  be  located  using  the 
hardcopy  printout  of  the  RTP  treatment  plan.  The  X  and  Y  values  must 
be  converted  into  units  of  millimeters  and  stored  in  global  variables 
on  the  disc  containing  the  program.  The  RTP  software  stores  the  X  and 
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Y  coordinates  in  a  form  that  requires  division  of  the  value  by  16  to 
convert  to  units  of  millimeters.  Each  external  contour  is  then  stored 
with  its  X  coordinates  stored  in  one  global  variable  and  its  Y 
coordinates  in  another. 

Sorting  routines  must  be  used  to  go  through  the  values  of  the  X 
and  Y  coordinates.  The  program  requires  the  maximum  X  and  Y  values 
and  the  minimum  X  value  for  each  of  the  external  contours  for  later 
operations.  The  maximum  Y  values  for  each  of  the  external  contours 
are  sorted  to  find  the  maximum  Y  value  out  of  all  the  contours. 

The  maximum  Y  value  is  used  to  shift  the  Y  coordinate  system  so 
that  Y  is  equal  to  zero  at  the  maximum  point  of  the  patient  surface 
and  all  other  y  values  are  less  than  zero.  This  is  necessary  for  the 
router  driver  routine  so  that  the  router  cuts  into  the  filter  mold 
where  there  are  tissue  deficits.  The  X  coordinate  system  is  shifted 
by  subtracting  the  value  for  the  X  coordinate  of  the  central  axis 
point  from  all  the  X  values.  This  sets  X  equal  to  zero  on  the  central 
axis  (see  Figure  3). 

The  data  points  of  each  external  contour  are  accessed  and  those 
points  which  fall  below  a  line  drawn  from  the  minimum  X  to  maximum  X 
points  are  deleted  (see  Figure  4).  The  remaining  contour  coordinates 
are  then  stored  in  ranked  order  from  the  maximum  X  value  down  to  the 
minimum  X  value.  This  enables  the  router  to  begin  at  one  end  of  the 
filter  mold  and  cut  its  way  to  the  other  end  for  each  CAT  scan  slice. 
The  router  can  not  follow  Indentations  along  the  X  axis  so  those  must 
be  deleted  prior  to  the  router  driver  routine  because  the  router  only 
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Figure  3.  The  figure  shows  the  coordinate  system  at  the  styrofoam 
filter  block.  The  dimensions  of  the  filter  are  10  inches  by  10 
Inches  by  3  Inches.  Point  A  Is  the  central  axis  point  at  the  top 
of  the  filter.  This  is  the  origin  of  the  coordinate  system  at  the 


filter 


Y 


XMIN 


Figure  4.  The  figure  gives  an  example  of  the  deletion  of 
indentations  along  the  X  axis  and  shows  the  cutoff  line  for  the 
points  to  be  deleted  on  the  bottom  of  the  CAT  scan  external 
contour.  Beginning  at  point  XMAX  and  moving  counter-clockwise, 
each  X  coordinate  value  is  compared  with  the  next  X  value  on  the 
contour.  If  the  following  X  value  is  not  less  than  the  preceeding 
value  it  is  set  equal  to  the  preceeding  value  and  rewritten  into 
the  file.  This  causes  indentations  similar  to  A  to  be  replaced 
with  a  straight  line  such  as  B.  This  enables  the  router  to  track 
properly  over  the  contour  in  the  X  and  Y  dimensions.  When  point 
XMIN  is  reached  the  coordinates  between  XMIN  and  XMAX  are  deleted 
along  with  their  corresponding  Y  coordinates. 
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moves  in  the  X  and  Y  dimensions  while  processing  a  CAT  scan  slice. 

Upon  completion,  it  then  increments  along  the  Z  axis  to  the  next  slice 
and  starts  over. 

The  X  coordinate  values  for  the  points  on  each  external  contour 
must  be  trigonometrically  reduced  for  the  reduction  in  lateral 
dimensions  at  the  filter.  The  diameter  of  the  source,  the  distance 
from  the  source  to  the  patient  at  the  central  axis  point,  the  distance 
from  the  source  to  the  filter  base  and  the  diameter  of  the  filter  base 
are  values  used  in  this  calculation.  The  source  is  then  treated  as  a 
virtual  point  source  located  at  increased  source-to-patient  and 
source-to-f liter  distances,  with  the  original  diameter  of  the  source 
acting  as  a  stop  (see  Figure  5). 

For  the  new  filter  holder: 

TAN-116 . 5/651-126 . 2/D2 

D2-126.5*(651/116.5)-707  (D2  is  updated  source-to-f liter 
base  distance) 

707-651-56,  80CH-56-856  (The  updated  source-to-patient 
distance) 

TAN-126. 5/707-R2/856 
R2-153mm 

This  makes  30cm  the  geometric  limitation  of  the  treatment  area 
along  the  X  axis  due  to  the  filter  holder  sides  interfering  with  the 
beam  divergence.  This  means  the  maximum  number  of  slices  that  can  be 
treated  with  the  new  filter  holder  in  place  is  twenty-nine.  The  size 
of  the  treatment  area  at  the  patient  can  be  adjusted  at  the  source 
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with  the  trimmer  bars  (which  limit  the  treatment  beam  divergence). 

This  allows  the  beam  divergence  to  be  limited  so  that  it  does  not 
contact  the  sides  of  the  filter  holder. 

Tangents  can  now  be  easily  calculated  for  transposing  the  X 
coordinates  at  the  patient  to  those  values  needed  at  the  filter. 

X(f liter )*(X( patient )/(856-Y(patient)))*651 

The  Z  coordinate  values  are  those  arising  from  the  one  centimeter 

thickness  of  each  CAT  scan  slice.  The  values  along  the  Z  axis  must 

also  be  reduced  to  those  needed  at  the  compensation  filter.  Since  the 

reduction  depends  on  the  distance  from  the  source  to  the  patient,  an 

average  distance  for  each  CAT  scan  slice  must  be  calculated.  This  is 

done  by  taking  one  half  of  the  height  of  the  CAT  scan  slice  along  the 

Y  axis  and  adding  this  to  the  source  to  patient  distance.  The  new 

source-to-patient  distance  is  used  in  calculating  the  tangent  values 

used  to  transpose  the  Z  coordinate  values.  The  one  centimeter  width 

of  each  slice  is  also  transposed  for  use  as  the  width  of  the  router 

bit.  The  Y  coordinate  values  are  transposed  into  the  values  needed  at 

the  filter  through  division  by  the  compensator  tissue  ratio.  This 

ratio  is  the  ratio  of  the  attenuations  of  the  compensator  and  the 

tissue  at  1.25  MeV,  the  average  energy  for  the  photons  from  the  Cobalt 

60  Therapy  Radiation  Source.  The  body  is  assumed  to  be  without  voids 
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and  with  a  uniform  density  of  lgm  /  cm  . 


The  router  driver  routine  is  set  to  start  with  transposed  X  and  Y 
coordinate  data  from  a  CAT  scan  slice  at  one  end  of  the  treatment 
area*  It  begins  at  the  maximum  X  position  for  the  slice.  The  router 
bit  width  is  given  and  the  bit  must  now  be  attached.  The  machine  then 
drives  the  router  bit  along  the  contour  given  by  the  X  and  Y  data 
points.  The  router  bit  is  then  removed  and  the  router  then  moves  to 
the  maximum  X  position  for  the  next  slice.  This  procedure  is  repeated 
until  the  filter  mold  is  completed. 

Since  there  is  no  router  machine  built  yet,  a  hardcopy  printout 
which  can  be  used  to  make  an  aluminum  block  filter  is  desired.  This 
filter  is  composed  of  half  inch  by  half  inch  blocks  of  aluminum  of 
different  thicknesses  stacked  to  the  proper  height  for  each  array 
position.  The  printout  is  ten  inches  by  ten  Inches  and  is  divided 
into  half  inch  by  half  inch  squares  along  the  X  and  Z  axes  (see 
Appendix  E).  The  Y  values  that  fall  within  each  of  the  squares  are 
averaged  and  their  mean  is  printed  within  the  square.  The  Y  value 
tells  one  how  tall  the  blocks  are  to  be  for  that  square. 


Figure  5.  The  figure  shows  how  the  source  will  be  treated  as  a 
virtual  point  source  and  the  limitation  of  the  maximum  treatment 
beam  divergence  due  to  the  filter  holder  itself.  Side  1  is  an 
example  of  the  Heustis  filter  holder.  The  angle  is  the  maximum 
half-angle  beam  divergence,  X2  is  the  half-width  of  the  filter,  XI 
is  the  half-width  of  the  source,  D1  is  the  length  to  be  added  to 
the  original  source-to-patient  and  source-to-f liter  distances  in 
order  to  treat  the  source  as  a  point  source  and  R1  is  the 
half-width  of  the  patient  that  can  be  treated  with  the  holder  in 
place.  Side  2  has  the  new  filter  holder  in  place.  D2,  the 
half-angle  divergence  and  the  half-width  at  the  patient  are 
different  from  the  corresponding  values  on  side  1  because  of  the 
different  original  source-to-f liter  distances. 
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III.  Validation 


The  router  driver  routine  cannot  be  validated  quantitatively.  The 
routine  does  drive  the  the  probe  on  the  vater  phantom  in  what  appears 
to  be  the  proper  manner.  The  X  and  Y  coordinates  when  printed  out  and 
the  Z  coordinates  and  the  width  at  the  filter  for  each  of  the  CAT  scan 
slices  are  correct  when  one  works  out  the  geometric  transformations 
done  to  them. 

Early  tests  of  the  program  used  external  contours  of  regular 
geometric  figures  of  different  sizes  (tilted  squares)  which  were 
manually  entered  using  existing  CMS  software.  The  program  output  for 
these  was  also  validated  geometrically. 

The  program  was  run  using  twenty-one  CAT  scan  slices  taken  of  the 
Rando  Phantom.  The  Rando  Phantom  is  a  human  skeleton  encased  in 
tissue  equivalent  material  (see  Figure  6).  The  hardcopy  printout  was 
used  to  build  an  aluminum  block  filter  (see  Figure  7). The  filter 
needed  a  special  holder  in  order  to  be  used  on  the  Cobalt  60  machine 
at  the  Medical  Center.  The  holder  was  made  out  of  aluminum  at  the 
AFIT  machine  shop  (see  Figure  8).  The  new  holder  has  both  a  rack  for 
the  shadow  tray  and  a  rack  for  the  filter.  The  AFIT  machine  shop  also 
made  the  half  inch  by  half  inch  aluminum  blocks  in  thicknesses  of  1, 

2,  5,  10,  20  and  30  millimeters.  The  aluminum  blocks  were  held 
together  using  rubber  cement  so  they  could  be  reused. 

The  phantom  was  then  treated  both  with  and  without  the  filter  in 
place.  Film  dosimetry  was  taken  during  each  treatment.  The  film 
dosimetry  was  then  evaluated  using  the  CMS  computer  driven 


Phantom  was  also  used  as  the  patient  in  the  film  dosimetry  testing 


the  Heustis  holder  and  a  tray  for  the  tissue  compensator.  The 


Figure  9.  The  isodensitometer  output  from  the  film  dosimetry 
results  for  treatment  without  the  compensation  filter  in  place. 


22 


IV.  Conclusions  and  recommendations 


The  computer  program  designs  a  compensation  filter  which  Is 
geometrically  accurate.  The  proper  reductions  in  the  lateral 
dimensions  do  take  place.  The  film  dosimetry  and  the  isodensitometer 
evaluation  would  give  a  more  accurate  picture  if  the  Rando  Phantom 
were  not  curved  on  the  underneath.  This  causes  the  dose  to  vary  in 
uniformity  at  the  film  because  some  portions  of  the  beam  pass  through 
less  tissue  before  exiting  the  phantom.  The  problem  could  be 
eliminated  by  using  thermoluminescent  dosimeters  placed  within  the 
phantom  itself  at  a  uniform  depth.  The  dosimeters  would  give  an 
accurate  picture  of  the  dose  distribution.  Time  constraints  prevented 
the  use  of  these  dosimeters.  After  further  verification  using  the 
TLD's,  the  routing  machine  should  be  constructed  and  evaluated. 

The  program  can  also  be  Improved;  There  are  three  areas  where 
improvements  might  be  made.  The  first  is  to  attempt  to  take  into 
account  the  varying  densities  of  the  patient.  Contours  of  areas  with 
a  different  density  from  the  surrounding  tissue  can  be  made  using  the 
RTP  software.  The  second  area  of  improvement  is  to  take  into  account 
the  scattering  due  to  the  filter  and  the  filter  holder.  The  third  is 
in  decreasing  the  calculation  time  of  the  program. 
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The  CMS  water  phantom  is  part  of  a  dosimetry  system  used  for 
plotting  isodose  contours.  It  is  a  remotely  controlled  device.  The 
water  phantom  is  capable  of  moving  its  probe  along  any  two  of  its 
three  dimensions.  The  driver  routine  in  the  program  will  drive  the 
Water  Phantom  probe  first  along  the  X  and  Z  dimensions  to  its  proper 
position  prior  to  following  the  transposed  contour  coordinates  of  the 
CAT  scan  slice.  The  Water  Phantom  will  the  move  along  the  X  and  Y 
dimensions  as  it  follows  the  contour  coordinates  of  the  scan  slice. 
This  procedure  is  repeated  until  all  the  CAT  scan  slices  have  been 
processed . 


phantom  probe  mount.  The  movement  of  the  water  phantom 


Figure  12.  The  figure  demonstrates  the  axes  of  movement  for  the 
Water  Phantom  probe.  These  X,  Y  and  Z  coordinates  will  be  the 
same  as  the  coordinate  system  at  the  filter  top  previously  shown 
in  Figure  3. 
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Figure  13.  The  Water  Phantom  50  pin  input  Connector  is  shown  with 
a  table  of  the  pin  functions. 
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Appendix  B:  Short  Guide  to  MUMPS 

Mathematical  Operators; 

+  Addition 
Subtraction 
*  Multiplication 
/  Division 
\  Integer  Division 
**  Exponentiation 

Mathematical  Operations 

There  is  no  hierarchy  of  operations.  Operations  are  performed 
strictly  from  left  to  right  unless  parentheses  are  used. 

Lines  with  mathematical  operations  must  be  preceeded  with  the  SET 
command.  One  SET  command  can  be  used  for  more  than  one  mathematical 
operation. 

S  X-(27**2)*4,  Y-X*37 

Arithmetic  Comparison  Operators: 

■  Equal  to  *■  Not  equal  to 

<  Less  than  '<  Not  less  than 

>  Greater  than  ’>  Not  greater  than 

Go  To  Commands 

Reset  the  pointer  so  that  the  line  label  immediately  following  the 


Go  To  command  is  the  next  line  executed  by  the  program 


G  LINE2  (Go  to  LINE2) 

Can  also  be  used  to  reset  the  pointer  so  that  the  next  line 
executed  is  in  another  program. 


G  LINE2  PR0GRAM2  (Go  to  LINE2  in  program  PR0GRAM2 ) 
There  are  also  conditional  Go  To  commands. 


G  LINE2 :X*5  (Go  to  LINE2  if  X  is  equal  to  2) 


If  Commands 

If  the  statement  is  true  then  the  rest  of  the  line  is  executed,  if 
the  statement  is  not  true  the  pointer  drops  to  the  next  line  without 
executing  the  other  commands  on  the  line. 


I  X-25  S  Y-K  G  LINE2  (If  X  equals  25  set  Y  equal  to  K  and  go  to 
LINE2) 


Variables 


Local  variables  can  be  up  to  15  alphanumerics< 


Arrays 


Arrays  are  subscripted  variables.  The  numeric  subscript  must  have 
a  value  between  0  and  8388607.  There  can  be  more  than  one  subscript 


for  a  variable.  String  subscripts  are  also  allowed. 
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Appendix  C:  Program  Documentation 


Program  LQAD1 

Purpose:  To  Initialize  values  needed  for  the  operation  of  the 
entire  linked  program. 

The  program  requests  the  operator  to  enter  the  X  and  Y  coordinates 
of  the  central  axis  point,  source-to-patlent  distance,  the 
source-to-f liter  distance,  the  number  of  the  CAT  scan  slice  containing 
the  central  axis  point,  the  compensator  tissue  ratio,  a  zero  If  the 
router  driver  routine  is  desired,  a  zero  If  the  filter  Is  on  the  top 
rack  and  the  patient  ID  number.  The  program  then  sets  N  equal  to 
zero.  The  program  then  increments  N  and  requests  the  contour 
description  for  the  CAT  scan  slice.  This  procedure  is  repeated  until 
all  of  the  descriptions  have  been  entered  and  assigned  an  N  value. 

The  operator  enters  a  zero  if  the  last  description  has  been  entered 
and  the  program  transfers  control  to  program  EXT11. 

Programs  EXT11  and  EXT22 

Purpose:  To  access  the  data  points  for  the  external  contours, 
convert  the  values  into  units  of  millimeters  and  store  them  In 
specific  global  arrays. 

Program  EXT11  sets  N  equal  to  zero.  The  value  for  N  is 
Incremented  and  the  program  accesses  the  MPF  to  retrieve  the  first 
block  of  the  file  containing  the  external  contour  stored  under  the 
contour  description  associated  with  N.  This  block  of  information  is 
placed  in  file  FID.  The  program  now  opens  file  FID3  connected  to 
global  ABC(N, "XCOMP" )  and  file  FID4  connected  to  global  ABC(N,”YCOMP") . 


INPUT  *  COORDINATE 


I  (PVT  T  COORDINATE 

~  ~~  ~  r  — 

I  (PUT  SOURCE -TO-PATIENT 
DISTANCE 

muz _ i: 

I MPUT  SOURCE-T  ‘-FILTER 
01  STANCE 

_  i" :  ~ 

INPUT  NUNSER  OP  SLICE 

containinc.  central  axis 

POINT 


INPUT  COKPENSATOR/TISSUE 
RATIO 

- - 1 - -  ' 

INPUT  0  FOR  ROOTER,  INPUT 

1  FOR  ALUMINUM  BLOC*  FILTER 

.  i  . 

INPUT  0  POR  FILTER  ON  TON 

RACK.  INPUT  1  POR  FILTER  ON 
BOTT*  RACK 


|  INP'-IT  pa-:ent  IP 


PLACE  FIRST  BLOC*  OP 

coordinate  values  in 

PILE  F ID 

READ  FIRST  X  VALUE 


CONVERT  X  COORDINATE  INTO 

millimeters  and  place  in 
PILE  PID3 

1 NCREMENT  TO  NEXT  X 


TES- 


STORE  FILE  FID 3  IN  GLOBA1 
ARRAY  ABC  N.'XCW) 


Figure  14. 


Linked  Program 


Flow  Chart  (Sheet  1 
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PUCE  2nd  BLOCK  OP 
COORDINATE  VALUES  III 

riu  pip 

BEAD  PIRST  X  VALUE 

- 

CONVERT  K  COORDINATE  IBTO 
MILLIMETERS  AMD  PLACE  III 
PILE  PID3 

INCREMENT  TQ  NEXT  K 


STCRE  PILE  PID3  IB  GLOBAL 
ARRAT  ABC(R."XCflM>-) 


PLACE  PIRST  BLOCK  OP 
COORDINATE  VALUES  IN 
PILE  PID 


Figure  14.  Linked 


PUCE  Zr.i  BLOCK  OP 
COORDINATE  VALUES  IN 

PILE  PID 

READ  PIRST  T  VALUE 

CONVERT  T  COORDINATE  INTO 
MILLIMETERS  AMD  PUCE  IN 
PILE  PID* 

INCREMENT  TO  NEXT  T 

^/ij 

istV 

TES 

STORE  PILE  PIDA  IN  GLOBAL 
ARRAT  AiCCN,*TCOKP') 

Program  Flow  Chart  (Sheet  2  of  7) 
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PLACE  CONTENTS  OP  GLOBAL 
ARRAY  A»C(N,’*C«P">  in 
PILE  FIDS 

~~ —  r  — 

FIND  MINIMUM  *  VALUE 
AND  PLACE  IT  IN  PILE 
FIDS 

INCREMENT  N 


YES 


PLACE  CONTENTS  OF  GLOBAL 
ARRAY  ABC(N.-YCOMP-)  IN 
FILE  FID* 

i  - 

FIND  MAXIMUM  T  VALUE 
AND  PLACE  IT  IN  FILE 
FIDS 

INCREMENT  N 


TES 


PLACE  FILE  FID5  IN 
CLOBAL  ARRAY  ABCCMAX') 


Figure  14.  Linked  Program  Flow  Chart  (Sheet  3  of  7) 


INCREMENT  N 


CLAST  »>-TES 


-(E) 


INCREMENT  W,  WD-0 


LAST  0>TIS 


~© 


PLACE  CONTENTS  OP  GLOBAL 
A* EAT  ABC(N,"ECCKP")  III 
PILE  PID3 


TAKE  THE  X  COORDINATES 
BETWEEN  HAXIHUH  X  AND 
MINIMUM  X  AND  PLACE  THEM 
IN  FILE  FIDS 


PLACE  THE  COWES  PONDING  T 
COMPONENTS  INTO  FILE  FIT'S 


PLACE  FILE  FIDS  IN  GLOBAL 
ABCCELOOT  > 

PUCE  PILE  FID6  IN  GLOBAL 
ARRAY  ABCCELIMT-) 


COPT  CLOBAL  ANRAT 
ABCCELDOt')  mo  CLOBAL 
A* RAT  ABC(N.-XCOKP-) 

COPT  CLOBAL  ARRAT 

abccelimv)  mo  global 

ARRAT  ABC<N,'TCOKP") 


PUCE  CONTENTS  OF  GLOBAL 
ARRAY  A*C(N,'XCCKP")  INTO 
PILE  PID3 


SET  X  EQUAL  TO  VALUE  CP 
WDth  WORD  OP  PILE  FID3 


/usrV 

Ot  COORDINATE) 

— yts  — 

~  1 

STORE  PILE  PID3  IN  CLOBAL 

no] 

ARRAY  AJC(N,  ~XC0MP“  ) 

1  " 

SET  T  EQUAL  TO  VALUE  OP 
WDth  WORD  OP  PILE  PID) 


WRITE  T  INTO  WDth  WORD 
OP  PILE  FID3 


Figure  14.  Linked  Program  Flow  Chart  (Sheet  4  of  7) 
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INCREMENT  n 


INITIALIZE  VALUES  NEEDED 
EON  TRANSPOSITION  Of  X  AND 
T  VALUES  TO  THOSE  TO  BE  USED 
AT  THE  FILTER 


INCREMENT  N 


PLACE  CONTENTS  Of  CLOBAL 
ARRAT  A BC(N,'XC<W)  INTO 
FILE  f  ID3 

~  :  i  : 

SUBTRACT  X  VALUE  OF 
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FILE  FID4 
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VDth  W«D  OF  FILE  FID3 
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Figure  14.  Linked  Program  Flow  Chart  (Sheet  5  of  7) 
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INITIALIZE  VALOIS  NEEDED 
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Figure  14.  Linked  Program  Flow  Chart  (Sheet  6  of  7; 
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Figure  14.  Linked  Program  Flow  Chart  (Sheet  7  of  7) 
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The  number  of  contour  points  for  the  external  contour  is  read  from 
file  FID  at  the  WDth  word  location.  R(N)  is  set  equal  to  the  number 
of  contour  points.  RK  is  set  equal  to  ((2*R(N))+7).  This  is  to  be 
used  as  a  limiting  value  since  it  will  give  the  location  of  the  last 
data  point  of  the  external  contour. 

The  X  coordinates  of  the  external  contour  are  accessed  first.  The 
accessed  data  points  are  divided  by  16  which  converts  the  units  into 
millimeters  and  the  results  are  written  into  file  FID3.  If  at  the  end 
of  the  first  block  of  data,  the  limiting  value  has  not  been  reached, 
control  is  transferred  to  program  EXT22  where  file  FID  is  closed  and 
the  global  array  of  the  external  contour  is  reaccessed  and  file  FID  is 
reopened  containing  the  second  block  of  the  global  array.  The 
limiting  value  RK  is  adjusted  to  provide  the  proper  limit  in  the 
second  block.  When  the  limiting  value  of  RK  is  reached,  whether  in 
the  first  or  second  block,  file  FID3  is  put  into  the  global  array 
ABC(N, "XCOMP" ).  Control  is  now  transferred  to  line  YRUN  in  program 
EXT11  and  the  Y  coordinates  are  processed  in  the  same  manner.  After 
they  have  all  been  processed  and  placed  in  global  array 
ABC(N,”YCOMP") ,  control  returns  to  line  XI  in  program  EXT11.  Here  the 
value  of  N  is  updated  for  the  next  external  contour.  The  data  points 
are  processed  and  stored  in  global  arrays  and  the  process  repeated 
until  until  N  reaches  the  limiting  value  of  Nl.  Control  then  passes 
to  the  first  line  of  program  XYMAX. 


Program  XYMAX1 


Purpose:  To  find  the  maximum  X  coordinate  and  the  maximum  Y 
coordinate  for  each  external  contour  and  place  the  values  in  global 
array  ABC ("MAX"). 

First  N  is  set  equal  to  1.  The  global  array  containing  the  X 
coordinates  of  the  first  external  contour  is  accessed  and  the  values 
are  placed  in  file  FID3.  The  value  at  word  0  is  read  and  X  is  set 
equal  to  that  value.  Y  is  set  equal  to  the  value  at  word  1.  The 
values  of  X  and  Y  are  compared.  If  X  is  greater  than  or  equal  to  Y, 
the  word  number  is  incremented  and  Y  is  set  equal  to  the  new  value.  X 
and  Y  are  compared  again.  The  process  continues  until  either  X  is 
less  than  Y  or  the  limiting  value  of  R(N)  is  reached.  If  X  is  less 
than  Y,  then  X  is  set  equal  to  the  value  of  Y,  and  the  word  number  is 
incremented.  When  the  word  number  reaches  the  limiting  value  of  R(N) 
the  value  of  X  is  the  maximum  X  coordinate  for  the  external  contour. 
The  global  array  ABC("MAX")  is  now  accessed  and  the  value  of  X  is  is 
placed  in  the  array.  N  is  then  incremented  and  the  X  coordinates  of 
the  next  external  contour  are  processed  in  the  same  manner.  After  the 
last  external  contour  has  had  its  X  maximum  determined,  control  the 
transfers  to  line  Y  where  N  is  reset  to  1  and  the  same  procedure  is 
repeated  to  find  the  maximum  Y  coordinate  values  for  the  external 
contours.  After  all  these  have  been  processed,  the  files  are  closed 
and  the  unnecessary  variables  are  killed.  Control  then  transfers  to 
line  one  of  program  XMIN1 
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Purpose:  To  find  the  minimum  X  coordinate  value  for  each  external 
contour  and  place  those  values  in  global  array  ABC("MAX"). 

N  is  set  equal  to  1.  The  global  array  containing  the  X 
coordinates  for  the  first  external  contour  is  accessed  and  the  values 
are  placed  in  file  FID3.  X  is  set  equal  to  the  value  at  word  0.  Y  is 
set  equal  to  the  value  at  word  1.  The  values  of  X  and  Y  are 
compared.  If  X  is  less  than  or  equal  to  Y,  the  word  number  is 
incremented  and  Y  is  set  equal  to  the  new  value.  X  and  Y  are  compared 
again.  If  X  is  still  less  than  or  equal  to  Y,  the  same  procedure  is 
followed  again.  If  X  is  gteater  than  Y,  then  X  is  set  equal  to  the 
value  of  Y,  the  word  number  is  incremented  and  Y  is  set  equal  to  the 
new  value.  This  procedure  continues  until  the  limiting  value  of  R(N) 
is  reached.  At  this  point,  the  value  of  X  is  the  value  of  the  minimum 
X  coordinate.  The  global  array  ABC(”MAX")  Is  accessed  and  the  value 
of  X  is  stored  within  the  array.  N  is  then  incremented  and  the  next 
contour  is  processed  the  same  way.  When  all  the  contours  have  been 
processed,  the  files  are  closed  and  the  unnecessary  variables  are 
killed.  Control  now  transfers  to  the  first  line  in  program  YMAX1. 

Program  YMAX1 

Purpose:  To  sort  through  the  previously  determined  maximum  Y 
coordinate  values  and  find  the  largest  maximum  Y  value.  This  value  is 
used  to  transpose  the  Y  coordinate  system  so  that  Y  coordinate  values 
are  equal  to  0  at  the  maximum  Y  point. 


The  global  array  ABCC’MAX")  is  accessed  and  its  value  placed  into 


file  FID5.  X  is  set  equal  to  the  value  of  the  maximum  Y  value  for  the 
first  contour.  The  word  is  incremented  and  Y  is  set  equal  to  that 
value.  If  X  is  greater  than  or  equal  to  Y,  then  the  next  word  value 
is  read  from  the  file  and  Y  is  reset  to  that  value.  X  and  Y  are 
recompared,  following  the  same  procedure,  unless  X  is  less  than  Y. 

When  X  is  less  than  Y,  X  is  reset  to  the  value  Y  and  Y  is  set  equal  to 
the  value  of  the  next  word  in  the  file.  When  the  value  of  the  last 
contour  has  been  processed,  YMAX  is  set  equal  to  X,  which  is  the  value 
of  the  largest  Y  coordinate  of  all  the  contours.  YMAX  is  then  stored 
in  global  array  ABC( "MAX" ) . 

The  global  array  ABC(N, "YCOMP" )  is  accessed  for  the  first 
contour.  The  values  are  placed  into  file  FID3.  Each  of  these  values 
is  read  and  has  the  value  YMAX  subtracted  from  it  and  then  is  replaced 
into  file  FID3.  After  all  the  Y  coordinates  have  been  processed,  file 
FID3  is  stored  in  global  array  ABC(N, "YCOMP" ) .  N  is  then 
incremented.  All  the  Y  coordinates  for  all  of  the  external  contours 
are  processed  the  same  way.  Then  all  files  are  closed  and  unnecessary 
variables  are  killed.  Control  is  transferred  to  the  first  line  of 
program  ELIM1. 

Program  ELIM1 

Purpose:  To  eliminate  the  external  contour  points  which  are 
underneath  the  patient  and  therefore  are  not  needed  for  the  filter 
design. 

The  program  accesses  the  global  array  ABC(N,"XCOMP")  for  the  first 
contour  and  places  it  on  file  FID3.  Next  the  program  accesses  the 
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global  array  ABCC'MAX”)  and  places  those  values  In  file  FID4.  The 
minimum  X  coordinate  value  Is  read  for  the  contour  being  processed  and 
X  is  set  equal  to  that  value.  The  program  then  reads  the  maximum  X 
coordinate  value  and  sets  Y  equal  to  it.  Z  is  set  equal  to  the  word  0 
value  of  the  stored  X  coordinates.  Z  is  then  compared  to  X  and  Y.  If 
Z  is  not  equal  to  either  one,  then  the  word  number  is  incremented  and 
Z  is  reset  to  the  new  value.  The  comparisons  are  made  again.  When  Z 
is  found  to  be  equal  to  Y  (the  maximum  X  coordinate),  MAXI  is  set 
equal  to  the  word  number  the  value  of  Z  was  taken  from.  The  word 
number  is  then  incremented  and  the  comparisons  continue  in  the  same 
manner  until  Z  is  found  to  equal  X  (the  minimum  X  coordinate),  MINI  is 
set  equal  to  the  value  of  the  word  number  Z  was  read  from.  The  word 
number  is  then  incremented  and  Z  is  then  recompared  with  X.  When  Z  is 
no  longer  equal  to  X  the  control  transfers  to  line  L3. 

The  global  array  ABC("ELIMX")  is  accessed  and  the  contents  placed 
in  file  FID5,  global  array  ABC("ELIMY")  is  accessed  and  the  contents 
placed  in  file  FID6.  The  contents  of  both  files  are  set  equal  to  0. 
The  global  array  ABC(N,"YC0MP”)  is  accessed  and  the  contents  placed  in 
file  FID7.  The  word  numbers  for  files  FID3  and  FID7  are  set  equal  to 
MAXI  and  the  word  numbers  for  files  FID5  and  FID6  are  set  equal  to  0. 
The  value  for  word  MAXI  is  read  from  FID3  and  FID7  and  written  into 
word  0  in  files  FID5  and  FID6  respectively.  The  word  numbers  are 
Incremented  and  the  process  repeated  until  the  word  number  for  files 
FID3  and  FID7  reaches  the  value  of  MINI. 

After  the  value  of  MINI  has  been  reached,  the  files  are  replaced 
into  the  global  arrays  they  originated  from.  The  contents  of  global 


array  ABC(”ELIMX")  are  copied  into  global  array  ABC (N,  "XCOMP ’’)  and  the 
contents  of  global  array  ABC(”ELIMY”)  are  copied  into  global  array 
ABC(N,"YCOMP'’).  The  entire  procedure  is  repeated  for  each  of  the 
external  contours.  After  the  last  contour  has  been  processed,  the 
files  are  closed  and  the  unnecessary  variables  are  killed.  Control 
now  transfers  to  the  first  line  in  program  INDENT. 


Program  INDENT1 

Purpose:  To  remove  indentations  of  each  external  contour  along  the 


X  coordinate  axis  so  the  router  can  properly  follow  the  contours. 


The  global  array  ABC(N, ’’XCOMP”  )is  accessed  for  the  first  contour 
(N»l)  and  the  contents  are  placed  into  file  FID3.  X  is  set  equal  to 
the  value  from  the  location  of  word  0.  The  word  number  is  incremented 
and  Y  is  set  equal  to  the  new  value.  If  Y  is  less  than  X,  the  word 
number  is  incremented  and  Y  is  reset  to  the  new  value.  X  and  Y  are 
compared  again.  This  procedure  continues  until  Y  is  greater  than  X, 
then  Y  is  set  equal  to  the  value  of  X  and  the  new  value  for  Y  is 
written  into  file  FID3  in  the  location  the  previous  Y  value  was  read 
from.  The  procedure  continues  until  the  last  word  has  been  read. 

File  FID3  is  then  placed  into  the  global  array  ABC(N, "XCOMP").  N  is 
incremented  and  the  next  contour  is  processed.  After  all  contours 
have  been  processed,  all  files  are  closed,  unnecessary  variables  are 
killed  and  control  transfers  to  the  first  line  in  program  XZER01. 
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Program  XZER01 

Purpose:  To  transpose  the  X  coordinate  system  so  that  the  origin 
is  at  the  central  axis  point. 

The  global  array  ABC(N, "XCOMP" )  is  accessed  for  the  first  contour 
and  the  contents  placed  in  file  FID3.  Each  of  the  X  coordinate  values 
is  read  and  has  the  value  of  the  X  component  of  the  central  axis  point 
subtracted  from  it.  The  modified  X  coordinate  values  are  then 
returned  to  file  FID3.  When  all  the  coordinate  values  have  been 
transposed,  the  file  is  placed  back  into  global  array  ABC (N, "XCOMP”) . 

This  procedure  is  repeated  for  all  of  the  external  contours.  When 
the  last  contour  has  been  processed,  the  files  are  closed,  unnecessary 
variables  are  killed,  and  control  transfers  to  the  first  line  in 
program  XYTRAN1. 

Program  XYTRAN1 

Purpose:  To  transpose  the  X  and  Y  coordinate  values  at  the  patient 
to  those  needed  at  the  filter. 

If  the  filter  is  going  to  be  on  the  bottom  rack,  the  source-to 
filter  distance  is  updated.  Using  the  half-width  of  the  source,  the 
half  width  of  the  filter,  the  source-to-patient  distance  and  the 
source-to-f ilter  distance  to  solve  for  the  tangent  of  the  maximum 
half-angle  beam  divergence.  This  value  is  used  to  determine  the 
distance  to  extend  the  source-to-patient  and  source-to-filter 


distances  in  order  to  treat  the  source  as  a  virtual  point  source. 


The  global  array  ABC(N, "XCCMP" )  is  accessed  for  the  first  contour 
and  the  values  are  placed  into  file  FID3.  The  global  array 
ABC(N, "YCOMP" )  is  accessed  and  the  values  are  placed  into  file  FIDA. 

X  is  set  equal  to  the  first  value  of  file  FID3.  Y  is  set  equal  to  the 
first  value  of  file  FIDA.  D2  is  set  equal  to  the  distance  (along  the 

Y  axis)  from  the  central  axis  point  to  the  point  being  processed.  D5 
is  the  set  equal  to  the  distance  (along  the  y  axis)  from  the  virtual 
source  to  the  the  point  on  the  contour  being  processed.  The  tangent 
of  the  half-angle  subtended  at  the  patient's  surface  in  the  X 
coordinate  plane  from  the  central  axis  to  the  point  on  the  contour  is 
calculated.  This  tangent  value  is  used  to  transpose  the  X  coordinate 
value  at  the  patient  to  the  value  to  be  used  at  the  filter.  If  the 
value  of  the  X  coordinate  at  the  filter  is  outside  the  volume  of  the 
filter,  then  “.he  contour  point  is  deleted  and  the  next  set  of  values 
for  X  and  Y  are  processed.  The  value  for  the  Y  coordinate  is 
transposed  to  that  used  at  the  filter  through  division  with  the 
tissue/compensator  ratio.  D7  is  set  equal  to  the  distance  (along  the 

Y  axis)  from  the  source  to  the  top  of  the  filter.  The  value  of  the  Y 
coordinate  is  now  compared  to  the  maximum  diameter  of  the  filter.  If 

Y  is  larger  than  the  filter  diameter,  then  the  point  is  deleted  and 
the  next  set  of  values  for  X  and  Y  are  read  and  evaluated  in  the  same 
way.  If  both  the  X  and  Y  coordinates  are  within  the  filter  volume 
they  are  written  into  files  F).J3  and  FIDA  respectively.  After  all  the 
contour  points  have  been  processed,  the  remaining  words  in  the  files 
are  set  equal  to  0. 


A7 


The  rest  of  the  external  contours  are  processed  In  the  sane 
manner.  After  the  last  one  is  done,  all  files  are  closed,  unnecessary 
variables  are  killed  and  control  transfers  to  the  first  line  of 
program  ZTRAN1. 

Program  ZTRAN1 

Purpose:  To  transpose  the  Z  coordinate  values,  where  Z  represents 
the  10  mm  thickness  of  each  CAT  scan  slice  along  the  Z  axis. 

The  global  array  ABC("MAX")  is  accessed  and  placed  into  file 
FID5.  XI  is  set  equal  to  the  maximum  Y  coordinate  value  (for  the 
first  contour),  which  is  read  from  the  file.  XI  is  divided  by  the 
compensator /tissue  ratio  to  transpose  the  value  to  the  present 
coordinate  system.  XI  is  divided  by  2  and  YMAX(N)  is  set  equal  to  the 
value.  YMAX(N)  is  used  to  calculate  the  width  of  the  router  bit  for 
each  slice.  YMAXl(N)  is  the  value  of  the  maximum  depth  (along  the  Y 
axis  for  each  slice)  divided  by  2.  D2  is  set  equal  to  the 
source-to-patient  distance  plus  the  YMAX(N),  this  allows  one  transpose 
the  Z  coordinates  trigonometrically  using  the  values  D2,  the  distance 
from  the  source  to  the  midpoint  of  the  patient  and  D1  the  distance 
from  the  source  to  the  midpoint  of  the  filter. 

The  transposed  Z  values  are  calculated  for  the  central  axis  slice 
first.  Calculate  the  values  for  the  tangents  using  distance  D2  and  a 
pair  of  Z  coordinate  values  set  equal  to  X  and  Y  (not  X  and  Y 
coordinates).  For  the  central  axis  slice  these  values  are  0  and  5. 

The  calculated  tangent  values  are  used  with  the  distance  D1  to 
calculate  the  transposed  Z  coordinate  values  XI  and  X2.  Twice  the 


48 


difference  between  X2  and  XI  is  the  value  assigned  to  WIDTH(N),  which 
is  the  width  determined  for  the  router  bit.  For  the  central  axis 
slice  the  Z  coordinate  location  for  the  center  of  the  router  bit  is 
Z(N)  equal  to  0. 

This  procedure  is  repeated  for  the  other  external  contours  with 
the  values  for  X  and  Y  (not  X  and  Y  coordinates)  incremented  by  ten 
after  each  slice  is  processed.  The  WIDTH(N)  value  is  the  difference 
between  X2  and  XI  (the  transposed  Z  coordinates).  The  value  of  Z(N) 
is  equal  to  the  X2  value  from  the  previous  slice  added  to  one  half  the 
value  of  the  difference  between  the  current  values  for  X2  and  XI. 

After  all  of  the  UIDTH(N)  and  Z(N)  values  have  been  calculated  for  all 
the  slices  the  unnecessary  variables  are  killed  and  the  file  is  closed. 

The  output  channel  to  the  printer  is  opened  and  the  proper  title 
is  printed.  The  values  for  all  the  WIDTH(N)  and  Z(N)  are  truncated  to 
3  decimal  points  and  printed.  Control  transfers  to  the  first  line  in 
program  WATER1. 

Program  WATER1 

Purpose:  To  drive  the  water  phantom  probe  in  the  same  manner  as 
the  tip  of  the  router  would  be  driven. 

If  the  value  of  A  is  equal  to  1,  this  means  the  aluminum  block 
filter  is  to  be  designed.  This  causes  control  to  transfer  to  the 
first  line  of  program  ALUM1.  If  A  is  not  equal  to  1,  program  WATER1 
is  followed. 
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The  limiting  variables  are  defined,  some  of  which  are  accessed 
from  global  IDV.  The  probe  of  the  water  phantom  is  positioned  at  the 
center  of  the  coordinate  system  on  the  filter  block  and  the  operator 
strikes  return  to  set  the  origin.  The  origin  is  the  central  axis 
slice,  at  the  central  axis  point  with  Y  equal  to  0. 

N  is  set  equal  to  1.  The  X  coordinates  of  the  first  slice  are 
accessed  and  placed  into  file  FID3,  the  Y  coordinates  are  accessed  and 
placed  into  file  FID4.  Z  is  set  equal  to  Z(N).  The  first  X  and  Y 
coordinates  are  read  and  the  probe  will  move  to  the  point  X,Y,Z(N). 

The  operator  must  strike  return  for  the  probe  to  be  driven  along  the 
external  contour  for  the  first  slice.  When  completed  the  operator 
must  strike  return  to  reset  the  probe  position  for  the  next  slice  to 
be  processed.  This  procedure  is  repeated  until  the  last  contour  has 
been  processed.  The  program  now  comes  to  a  stop. 

Programs  ALUM1,  HARD1  and  S0RTX1 

Purpose:  To  generate  a  hardcopy  design  for  the  aluminum  block 
filter  in  the  form  of  a  two  dimensional  array.  The  array  is  a  10  inch 
by  10  inch  square  divided  into  1/2  by  1/2  inch  squares.  Each  of  the 
1/2  inch  by  1/2  inch  squares  is  to  have  printed  within  it  the  height 
in  millimeters  of  the  compensating  material. 

The  Z(N)  values  are  converted  into  units  of  half-inches.  R1  is 
set  equal  to  the  distance  between  the  first  slice  and  the  central  axis 
slice,  R2  is  set  equal  to  the  distance  between  the  central  axis  slice 
and  the  last  slice.  R1  and  R2  are  raised  to  the  next  Integer  value. 

R1  is  set  equal  10  minus  Rl.  R1  is  now  equal  to  the  number  of  rows 


which  need  no  compensation  material  (up  to  the  midpoint  of  the 
filter).  The  20  R(K)  values  are  set  equal  to  0.  Control  is 
transferred  to  the  first  line  in  program  HARD1. 

The  output  port  to  the  printer  is  opened  and  the  title  and  the 
patient  ID  are  printed  out.  The  top  asterisk  boundary  is  printed  and 
the  Y(K)  values  for  each  1/2  inch  by  1/2  inch  square  are  printed. 
Control  transfers  back  to  line  L5  in  program  ALUM1.  The  value  of  R  is 
Incremented  and  control  then  transfers  back  to  line  LI  in  program 
HARD1,  where  the  array  row  of  the  filter  is  printed  out  with  its  Y(K) 
equal  to  0  values.  This  procedure  continues  until  R  reaches  the 
limiting  value  of  R1  plus  1.  After  all  the  empty  rows  in  the  first 
half  of  the  filter  have  been  output,  the  calculation  of  the  array  rows 
with  non-zero  Y(K)  values  can  begin. 

Z1  is  set  equal  to  the  value  of  R1  minus  10,  Z2  is  set  equal  to  Z1 
plus  1  and  Z3  is  set  equal  to  R2.  A  sorting  routine  is  begun  to 
determine  which  Z(N)  values  are  between  Z1  and  Z2.  NO  is  the  first 
slice  with  a  Z(N)  value  between  Z1  and  Z2,  NF  is  the  last  slice 
between  Z1  and  Z2.  Z1  and  Z2  are  now  incremented  and  control  now 
transfers  to  the  first  line  in  program  S0RTX1. 

The  global  arrays  ABC(N,"XC0MP")  and  ABC(N, "YCOMP" )  are  accessed 
and  placed  into  files  FID3  and  FID4  (for  contour  NO).  The  values  in 
FID3  are  converted  into  units  of  half-inches  through  division  by 
12.5.  A  sorting  routine  is  begun  to  determine  which  X  coordinates  lie 
within  the  boundaries  of  each  of  the  half-inch  squares  on  the  array 


Beginning  with  XI  set  equal  to  9  and  X2  set  equal  to  10,  the  first 
X  coordinate  is  read  (  the  maximum  X  coordinate).  If  X  is  greater 
than  X2  the  point  is  deleted  and  the  next  point  is  read.  If  X  is  less 
than  X2  but  not  between  XI  and  X2,  then  both  XI  and  X2  are  decremented 
by  1.  X  is  recompared  to  the  new  values.  When  X  falls  within  the 
boundaries,  P0  is  set  equal  to  the  word  number  the  X  coordinate  was 
read  from.  This  continues  until  X  no  longer  lies  within  the 
boundaries  of  XI  and  X2.  PF  is  set  equal  to  the  word  number  this  X 
coordinate  was  read  from.  Now  the  corresponding  Y  coordinate  values 
from  locations  P0  through  PF  have  their  values  multiplied  by  a  -1  and 
are  averaged  together,  with  Y(N,K)  set  equal  to  the  mean.  The  N  is 
the  number  of  the  slice  the  values  came  from  and  K  is  the  number  of 
the  square  the  values  fall  within.  Some  of  the  X  coordinates  make 
large  jumps  in  value  so  that  they  do  not  fall  within  the  boundaries  of 
one  or  more  squares.  These  jumps  occur  where  there  are  flat  spots  on 
the  contour  (where  the  Y  coordinate  values  don't  change).  This  would 
cause  false  values  of  0  to  be  assigned  to  Y(N,K).  The  program 
compensates  for  this  by  assigning  the  previous  Y(N,K)  value  to  the 
squares  effected. 

This  process  continues  until  all  the  Y  components  have  been 
assigned  to  Y(N,K)  variables.  This  is  done  for  CAT  scan  slices  NO 
through  NF.  The  Y(N,K)  values  for  each  K  value  are  averaged  and  Y(K) 
is  set  equal  to  the  result.  Control  then  transfers  to  line  LI  in 
program  HARD1.  The  asterisk  boundary  is  printed  and  the  Y(K)  values 
are  printed  in  the  array  row.  Control  then  transfers  back  to  line  L5 
in  program  ALUM1. 


The  next  values  for  variables  NO  and  NF  are  calculated  In  the  same 
manner  as  before.  Control  transfers  to  program  S0RTX1  where  the  next 
set  of  Y(K)  values  are  calculated.  Control  transfers  to  program  HARD1 
where  the  new  row  Is  printed  out.  This  continues  until  all  the 
contour  slices  have  been  processed.  The  Y(K)  values  are  then  reset 
equal  to  0  in  program  ALUM1  and  control  transfers  back  to  program 
HARD1  where  the  rest  of  the  array  rows  are  printed  out.  The  program 
then  stops. 


Appendix  D:  The  Linked  Program 


LOAD  I '  00 'SO  I 


01  DEC  93  10  05  10  RTP,  JUS 


.PROGRAM  FOP  OPERATOR  EN'ERED  Data 

P  "ENTER  VAuUE  IN  mm  S  POP  X  COMPONENT  OP  CEn'Ral  ai 15  PT  * .  '.xi  1 

P  "ENTER  VAlUE  IN  MM  S  POP  V  COMPONENT  OP  CENTRAL  ACS  PT',  '.VI.  ' 

P  "ENTER  DISTANCE  PROM  SOURCE  TO  CENTRAL  AXIS  PT".  '.Dl,  ' 

R  "ENTER  DISTANCE  FROM  SOURCE  TO  BASE  OF  TOP  FILTER".  '.03  1 

R  "ENTER  TmE  number  OF  TH£  CAT  SCAN  SLICE  CONTAINING  THE  central  AXIS  PT", 
P  "EN'EP  THE  COMPENSATOR /T ISSUE  RAT  10" .  ' . CTR,  ' 

P  "ENTER  0  PGP  ROUTER.  ENTER  I  POP  ALUMINUM  BLOCK  FILTER".  '  .  A.  ' 

P  "ENTER  0  IF  FILTER  ON  TOP  RACK,  l  IF  FILTER  ON  BO'TOM  RACK".  ’ ,  T,  • 

P  "ENTER  PATIENT  ID".  '.PID.  • 

S  N.Q 
S  N-N-l 

R  ENTER  CONTOUR  OISCPIPTION,  IF  LAST  SLICE  HAS  BEEN  ENTERED. ENTER  0".  '  .  CDI 
I  CDES < N  > »0  S  NI*N  S  OUIT 
G  LI 

0  X'EXTIt 
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•'.3  -PROGRAM  TO  PINO  XMAX  &  YMAX  IN  EACH  EXTERNAL  CONTOUR  AND  PLACE  THEM  IN  FIDS 

ZC  * 

S  N«1 

XCOMP  S  I«-A8C(N.  "XCOMP") 

ZO  "XCOMP".  FIDO 
ZC  FI 03  1 
S  .-0 
S  L»1 

ZR  FI 03  «J  X 
-AX  ZR  FI03  *L  Y 

C  A  X '<Y  G  8 

»  S  L«L*1  G  C  L«R:n)  C  -AX 

3  S  X»V  S  L“L*1  3  C  L«R  •<>  0  MAX 

C  G  ST  N»1  C  W 

ST  S  I»'ABC( "MAX")  Z A  "MAX"  1  CO 

ZO  “MAX"  FIDS 
ZC  FI09  1 

F  1-0  1  299  ZU  FIDS  I  0 

s  wo-a 

U  ZU  FIDS  UD  X 

S  N-N*l ,  UD-UD+1 
C  V  N-Nl  ZC  FIDO  0  XCOMP 
>  S  N«1 

•CC«F  S  l-"A3C<N,  "YCOMP") 

ZO  "YCOHP"  FID4 
ZO  FI 04  l 
'  S  J-0 
S  L«1 

ZR  FI 04  *J  X 

-AX2  ZR  FI 04  #L  Y 

C  A2  X>Y  G  82 

S  L-L*l  0  C2  L-R(N)  G  MAX 2 
S  X-Y  S  L-L-l  C  C2  L»R(N)  G  MAX2 
G  ST2  N-l  C  U2 
2  S  UD-N1-1 

ZU  F109  MO  X 
S  N-N*l.  UO-UD-l 
G  0  N-Nl  ZC  FI 04  C  YCOMP 
3  S  I-'A»C<"HAX") 

ZF  FI09  I 

3UIT  ZC  •  K  X.Y.  I.FI03.  J.  L.  FIDS,  UD.F1D4  9  X4~XMJNl 


Figure  15.  Linked  Progams  (Sheet  4  of  14) 
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'.;7  ;  PROGRAM  TO  ELIMINATE  THOSE  PT5  BELOW  THE  XMAX  \  x"IN  PT3  UNDER  THE  PATIENT 

ZC  •  3  N«0 

s  wo;«2»n;-2 

.1  S  N»N*’..Q»0  S  QUIT  N-Nl 

S  r«*ABC(N,  "XCOMP  • :  ZA  "XCOMP"  1  CO 
ZO  "XCOMP-  F 103 
ZC  F103  l 

S  I«*A3C  <  "MAX  ’  )  ZA  -'MAX'-  1  CO 

ZO  -MAX-.FI04 

ZC  F 104  1 

I  N> 1  S  UOl-UDl-l 

ZR  FI 04  *WOl  X 

S  UD2-N-1 

ZR  FID4,  -U02  V 

S  WD»0 

ZC  FI 04 

.2  ZR  F 1 03  -WD  Z 

I  Q-l  C  S2  Z-X  C  L3 

C  31  Z-Y  C  32  Z-X  S  WD-WD-1  C  L3  UD-*R  (N>  3  L2 

31  S  MAX  1 -WD  W0-W0*1  0  L3  WD-R(N)  C  L2 

32  S  MINi-WD,  Q-l.  WD-WD-1  0  L2 

-3  S  I-'A3C<"ELIMX")  ZA  “ELIMX"  1  CO 
ZO  "ELIMX"  FIDS 
ZO  FI  OS:  I 

S  I  ABC  (  "  EL  I  MY  -  >  ZA  "ELIMY"  1  CO 
ZO  "ELIMY"  FID6 
'  ZO  FI06  1 

F  1-0  1  2SS  ZW  FIDS:  10 
F  1-0  1  233  ZW  F 1 06  I  0 
S  I-‘A3C(N,  "YCOMP  )  ZA  " YCOMF " .  I  CO 
ZO  "YCOMF"  FID7 
ZO  FI 07  1 
S  WD— 0.  MAX -MAX  1 
.4  ZR  FID3  -MAX  X 

ZR  FI07  .MAX  V 
ZW  FIDS  WD  X 

ZW  FI 06  WO  Y  0  LS  MAX-MINI 
S  WD-WO-1.  MAX-MAX*!  0  L4 

-.3  S  WD < N ) -WO  ZF  FIDS  1.FID6  1 ,  F 103  1 .  FID7  l  ZC  FIDS.  FID6.  FID3.  FID7 
3  I-'ABCt "ELIMX")  ZA  "ELIMX". 1  CO 
ZO  "ELIMX"  FIDS 

S  I-'ABC<N.  “XCOMF")  ZA  "XCOMF"  1  CO 
ZO  "XCOMF"  FID3 
ZX  FIDS:  FIDO 
ZC  FI03.  FIDS 

8  I-^ABCI "ELIMY")  ZA  "ELIMY".  1  CO 
ZO  "ELIMY"  FID6 

3  I-^ABCIN. -YCOMF")  ZA  "YCOMF".  1:  CO 
ZO  “YCOMF":  FID4 
ZX  FI 06  FI 04 
ZC  FI 04.  FI06  0  LI 

WIT  K  I.R.FI03.FI04.  0.W02.  X.Y,  Z. MAXI. MINI, FIDS.  FID6.FI07.  MAX  0  XE^INOENTl 
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v.8  .PROGRAM  TO  REMOVE  INDENTATIONS  IN  THE  *  DIRECTION  PR  ON  The  OA'a  PTS 

ZC  •  S  N-»0 

LI  S  N»N*1.WD»0  £  SUIT  N«N1 

S  I«*A8C  :N.  "XC2MP  •' )  ZA  "XCOMP".  1  CO 
ZO  " XCOMP "  FIC3 
ZC  FI03  1 

L2  ZR  FID3  *UD  X 

S  UD-UD-I  0  L3  UD'UDIN) 

ZR  FI 03  »UD  Y 
C  L2  VO!  C  A1 

A 1  S  »*I  ZU  FI 03  HO  V  6  L2 
-3  ZP  FIDS  1  ZC  F103  6  LI 

QUIT  K  I.F103.  X.V  0  %9-XZEROl 
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\»  .  PROGRAM  TO  TRANSPOSE  X  COMPONENTS  SO  THAT  X-0  AT  CENTRAL.  AXIS  POINT 

ZC  *  S  N-0 

L.I  S  N»N*1  C  QUIT  N-Nl 

S  I»'ABC(N,  "XCOMP")  ZA  "XCOMP"  1  CO 
ZO  "XCOMP"  FIDO 
ZC  FI 03  1 
S  UD-0 

.2  ZR  FID3  »UD  X 

S  X«X-X1 
ZU  FI 03  UD  X 

S  UD« ! UD* 1 )  C  L3  UD»(U0<NM-1  >  Q  L2 
l.3  ZP  FID3  1  ZC  •  0  LI 
QUIT  ZC  *  K.  I.  FIDS.  X,  XI  C  %10~XYTRAN1 


Figure  15.  Linked  Progams  (Sheet  8  of  14) 
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•'.11  ,  PROGRAM  TO  TRANSPOSE  Z  COMPONENTS  .  Z  REPRESENTS  THE  10  MM  THICKNESS  OF  EACH  CT  SCAN 

zc  • 

;.l  S  I»’A8C<  MAX")  ZA  "MAX’  1  CO 

ZO  "MAX"  FIDS 
ZC  FIDS  1 
S  WD—N 1-1 
S  N-l 

L2  zr  FIDS  *UD  XI 

S  YMAX iN>*< X 1-VMAX > /CTR/2  .UPDATE  YMAX  t  N )  TO  PRESENT  COORDS 
S  YMAX1(N)»(75t(X1-YMAX)  IKC'R/a 
S  N-N-l 

S  WD«WD*1  C  L3  WD-< (2*N1 )-2>  5  L2 
-3  S  N=N2  .CENTRAL  AXIS  SLICE 

.4  S  0  1  =07-YMAX  >'  N  )  .  D2-03-.YMAX  1  (  N  > 

S  X  =0.  Y-S 

S  TAN1«X/D2  TAN2-Y/D2 

S  X 1«TAN1»DI .  X2-TAN2«D1.  X2<N)-X2,  Z(N)»0 
K  YMAX ;N>,  YMAX 1 (N> 

S  WIDTH<N)«(X2-X1 )»2.  N»N*1 
S  X-5.  Y«1S 

Li>  S  01*D7-YMAX:N>  02-DS*YMAXl<N) 

S  W-N-l 

S  TAN1-X/D2.  TAN2-Y/D2 
S  X1-TAN;»D1,  X2«TAN2*D1.  X2<N)«X2 
K  YMAX (N)  VHAX 1 (N) 

S  UIDTH(N)«X2-X1.  Z(N)«X2(U)*< (X2-X11/2).  N«N*1.  X-X  +  10.  Y»Y*10  C  L7  N-Nl  C  L6 
S  N-N2,  X-S.  Y-15 
S  N-N-l  C  QUIT  N-0 
S  W-N-l 

S  D1-07-YMAX.N).  02-D3-YNAX1 (N> 

S  TAN1-X/D2- TAN2-Y/D2 
S  Xl-TANl-Dl.  X2-TAN2-D1.  X2<N)-X2 
K  YMAX(N).  YMAX  1  < N  ) 

S  WIDTH(N)«X2-X1,  Z<N>  —  <X2<W>*( <X2-X1 >/2> >.  X-X-10.  Y-Y-10  0  LB 
QUIT  M  X.  XI,  X2.  Y.  Dl.  D2.  D3.  D7  W.  WO.  YMAX.  YHAXfN),  YHAX1  (N>,  TAN,  TAN1 .  TAN2  ZC  * 

0  2  PRT. 0  0 
S  N-0 

W  2.  "OUTPUT  FROM  PROGRAM  ZTRAN”,  !.  1 
L9  S  N-N-l  3  L10  N-Nl 

S  WIDTH(N)-WIDTH(N)»1000\1  .  Z  (  N  >  -  Z  (  N  )  *  1 000 \ 1 
S  WIDTH<N)»WIDTM<N> / 1000. Z<N1«Z(N)/1000 
W  :  2.  “N-".  N.  I>3.  "WIDTH-" .  WIDTH! N  > .  L3.  “Z*".  Z(N>.  ! 

I  A-l  K  WIDTH(N) 

C  u9 

L 10  K  I, FIDS, X2(N)  3  X12~WATER1 


Figure  15.  Linked  Progams  (Sheet  10  of  14) 
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i  PROGRAM  TO  OUTP'jT  DESIGN  FOR  ALUMINUM  SLOCK  FILTER 
S  N»C.  K«20.  Zl-O,  Z2-0.  Z3-1.  NO-O.  NF-0 
S  N«N*1 

S  Z  <N ) -Z  ( N)  / 12  S  ,  Z<N>  IN  .MTS  OF  HALF-INCHES 
C  L2  N-(Nl-l)  G  Ll 
S  R1-Z(N2)-Z< 1) .  R2«Z<N2)*Z<N) 

S  R1«<R1\1>*1.  R2-(R2\1)*1 
S  Rl-lO-Rl 
S  K-0,  R-l 
S  K-K-M 

S  Y(K)»0  C  L*  K-20  G  L3 
C  %19'HARDl 

C  L7  K-0  S  R-R*  1  G  LS  R-<R1»1>  G  Ll~HAR01 
S  N-0.  Zl-Rl-10.  Z2-Z  1*1.  Z3-R2.  U-0,  H-0 
S  N-N-l  G  LlO  Z1  <Z3  1  M»0  G  LlO  N-Nl 

I  Z<N)'<Z1  S  U-U-H  I  Z<N)'>Z2  S  M-H-l  0  FIRST  (1»|  I  N-(Nl-I)  S  NF-N  G  OUT 
0  LAST 

I  (1-1  S  NO-N.  W«1  G  L9  N-(Nl-l)  0  NEW 
I  W>M  S  NF-N-1  G  OUT 
0  L7 

X  ft*  1  S  NF>N 

S  W-0, M-0.  Zl-Zl-l.  Z2-Z2+1  0  X16AS0RTX1 
S  K-0 
S  K-K-l 

S  V ( K ) -0  G  L12  K-20  G  Lll 
G  L 1 "HARD1 
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Figure  15.  Linked  Progams  (Sheet  14  of  14) 
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The  Printout  (reduced)  for  the  Aluminum  Block  Filter 
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